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Abstract. Lattice quantum chromodynamics (QCD), defined on a discrete space time lattice, leads 
to a spectacular non-perturbative prediction of a new state of matter, called quark-gluon plasma 
(QGP), at sufficiently high temperatures or equivalently large energy densities. The experimental 
programs of CERN, Geneva and BNL, New York of relativistic heavy ion collisions are expected to 
produce such energy densities, thereby providing us a chance to test the above prediction. After a 
brief introduction of the necessary theoretical concepts, I will present a critical review of the exper- 
imental results already obtained by the various experiments in order to examine whether QGP has 
already been observed by them. 



1. Introduction 

As is well known [1], the standard model of particle physics, SU(3) C x SU(2) W x U(l)y 
broken spontaneously to SU(3) C x U(l) em , has been tested with great precision at LEP. 
All these tests rely heavily on the fact that the corresponding coupling is weak and hence 
the usual weak coupling perturbation theory can be employed in deriving the required 
theoretical predictions. Since the electromagnetic and weak couplings are indeed rather 
small in the currently accessible energy range, a em — 7.3 x 10~ 3 and a„. ~ 3.4 x 10~ 2 , use 
of perturbation theory is not a serious limitation in the precision tests of the electroweak 
theory. However, the strong interaction coupling, a s , is (i) a strongly varying function of 
energy in the same range, (ii) about 0. 1 1 at the highest energy at which it has been measured 
so far, and (iii) ~ 1 at a typical hadronic scale. Therefore, testing the strongly interacting 
sector of the standard model using only perturbation theory is a major shortcoming of the 
precision tests of the standard model. 

Formulating quantum chromodynamics (QCD), which is an St/ (3) gauge theory of 
quarks and gluons, on a discrete (Euclidean) space-time lattice, as proposed by Wilson 
[2], and simulating it numerically as first shown by Creutz [3], one can obtain [4] several 
post-dictions of QCD in the non-perturbative domain of large oc v . These include qualitative 
aspects, such as quark confinement and chiral symmetry breaking, and quantitative details 
such as hadron masses and their decay constants. While these agree with the known exper- 
imental results within the sizeable theoretical errors, it is fair to say that no serious exper- 
imental test of any non-perturbative prediction of QCD has so far been made. Relativistic 
heavy ion collisions offer a great window of opportunity to do so. Application of lattice 
techniques to finite temperature QCD has resulted in the prediction [5] of a new state of 
matter, called Quark-Gluon Plasma(QGP), at sufficiently high temperatures or energy den- 
sities. Chiral symmetry, broken spontaneously at zero temperature, seems to be restored in 
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this new phase characterised by a much larger degrees of freedom characteristic of "free" 
quarks and gluons. Nevertheless, the phase appears to be inherently non-perturbative in the 
experimentally interesting range of 1 < T/T c < 4-10, where T c ^150 MeV is the transition 
temperature at which the energy density varies most rapidly. The energy density, e, in this 
range is 15-20% smaller [6] than the value of the corresponding ideal gas of quarks and 
gluons whereas a maximum of 3-5% deviation is allowed for a weakly interactive perturba- 
tive QGP. While the precise values for e, or T c , as well as the nature of the phase transition 
(whether first order or second) depend on the number of light quark flavours, the quoted 
values above being for 2 flavours of mass about 15 MeV, many simulations with varying 
numbers of light flavours suggest that an energy density greater than 1 GeV/fm 3 is needed 
to reach the QGP phase. 

Collisions of heavy ions at very high energies can potentially produce regions with such 
large energy densities. Furthermore, since the transverse size of such regions is given by the 
diameter of the colliding nuclei, one can hope that these collisions will satisfy the necessary 
thermodynamical criteria of large volume (L <~ 2Ra » ^gcp) an d many produced particles. 
A crucial unanswered question is whether thermal equilibrium will be reached in these 
collisions, and if yes, when it will be reached and how. Many different attempts have 
been made, and are being made, to address these issues. Here we will follow Bjorken's 
picture as it is most widely used in the field. Bjorken argued [7] that for sufficiently high 
energies, yfs > 15A GeV where yfs is the total CMS energy of the two colliding nuclei 
of mass number A, the nuclei bore through each other and leave behind a baryonless blob 
of produced particles in the center (around y cm = \ In \(E + P L ) j(E — P L )] <~ 0). After an 
equilibration time To, the energy density in the blob was estimated by Bjorken to be 

J4 To ay 

where the effective area A — %R\ = 3.94 A 2 / 3 fm 2 and dEj/dy is the measured transverse 
energy per unit rapidity round y cm m 0.0. 

The Bjorken scenario for how the (thermally) equilibrated blob evolves is also the back- 
bone of the analyses seeking to extract information from the data on whether QGP did form 
in the heavy ion collisions. According to this scenario, the hot blob cools by expanding 
and the matter in it goes through various stages such as QGP, mixed phase and a hadron 
gas, depending on the initial energy density reached and the equation of state. A further 
rapid expansion of the hadron gas leads to such large mean free paths for the hadrons that 
they essentially decouple from each other. If this freeze-out is sufficiently fast, the free- 
streaming hadrons, %,k, ■ ■ ■ etc. will, however, retain the memory of the thermal state from 
which they were born by having thermal momentum distributions. Thus the information 
from observables related to light hadrons can tell us about the temperature at this 'thermal 
freeze-out' and the velocity of expansion. To get a glimpse at still earlier times, one has 
to turn to 'harder' probes which typically involve a larger scale such as masses of heavy 
quarkonia, as we will see below. 

The Bjorken scenario is for very high collision energies, when one expects to obtain a 
baryon-free region, and Eq. 1 is a valid description. The present collision energies may not 
be sufficient for it to hold, i.e., there may be a lot of baryons deposited in the central region 
of 

ycm ~ 0.0. A reliable analogue of Eq. 1 is however not available in that case. Note 
that even the theoretical estimate from lattice QCD above was for a baryonless case. In 
addition to temperature, one can also imagine increasing the baryon density of the strongly 
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interacting matter or equivalently increasing the baryonic chemical potential /j# and obtain 
a baryon-rich plasma. In principle, one knows how to handle the case of a nonzero baryon 
density on the lattice but it has so far turned out to be difficult in practice. Usual lattice 
techniques fail for nonzero /lib due to technical reasons [5] and attempts to overcome [4] 
these have not been successful either. No reliable lattice estimates are therefore available 
in that case. Using models based on underlying symmetries, it has been recently argued 
[8] that the T - /lib phase diagram of QCD with realistic quark mass spectrum should have 
a critical point at a nonzero /jb- The analysis of heavy-ion data by Ref. [8] did not reveal 
any such critical point. By varying the energy of the colliding beam of heavy ions, one 
may hope to unearth such a critical point. While upcoming experimental runs at CERN 
at y/s <~ 9 A GeV will look for such a critical point, the model considerations above are 
inadequate to provide reliable information on the energy density of the QGP phase for 
nonzero /jg. Thus a greater theoretical effort is required to firm up the QCD prediction for 
the energy density for nonzero /jg and also to obtain the analogue of Eq. 1 . Of course, one 
can in stead go for higher energies to test QCD, where one expects to obtain a baryon-free 
region, making both the lattice estimate and Eq. 1 more accurate descriptions. 

2. Results from CERN 

The experimental programs of high energy heavy ion collisions are being pursued actively 
at present in Brookhaven National Laboratory (BNL), New York and CERN, the European 
Laboratory for Particle Physics, Geneva. Au-Au collisions at s/s = 4.7 A GeV ~ 0.92 
TeV have been studied at BNL while Pb-Pb collisions at y/s = 17.3 A GeV ~ 3.6 TeV 
have been investigated at SPS, CERN using beams of gold ions at 2.1 TeV/c and lead ions 
at 32.9 TeV/c respectively. Earlier sulphur beam at 6.4 TeV/c was used on sulphur and 
uranium targets at SPS, CERN and those results form a benchmark over which several 
aspects of Pb-Pb collisions have been compared. I will focus largely on the latter since 
they correspond to the highest y/s used so far. Due to space restrictions, I will also have to 
restrict myself to highlights and I have to refer the reader for more details to the proceedings 
of Quark Matter conferences [9]. 

2.1 Initial Energy Density 

The NA49 experiment reported measurements on dEj jdx\ quite a while ago [10] and re- 
ported dEj /dy~ 405 GeV for Pb-Pb. Using a canonical guess of 1 fm for the formation 
time, one obtains from Eq. 1 

e™- pfe (lfm) = 2.94±0.3GeV/fm 3 , (2) 

which is certainly above the characteristic QGP-phase values from lattice QCD mentioned 
in sec. 1. Since appreciable numbers of baryons at y cm <~ have been observed at SPS, 
it is doubtful that the current energies are high enough for creating a baryon-free region 
assumed for Eq. 1. Nevertheless, Eq. 1 has been frequently used to estimate the energy 
densities achieved up to now. One has to be cautious therefore and make sure that other 
independent estimates are also similar and they do appear to be so. 
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2.2 Hadron Yields 

Assuming that a thermal freeze-out is triggered by a rapid expansion, one expects the 
momentum spectra of various hadrons to reflect the freeze-out temperature, 7> , which will 
be blue-shifted by the collective expansion. For small transverse momenta, pj <C m, we 
expect the inverse slope oftransverse mass distribution, da/dimr — m) withm^ = pj+m 2 , 
to be given by r s i ope = Tf + ^m{vj) 2 . Thus one expects, r s i ope to vary linearly with the 
mass of the observed particle, m. Fig. 1 displays the dependence of r s i ope on m for various 
particles produced in the Pb-Pb collisions in the different CERN experiments. A clear 
linear rise is evident for most of them except the heavier strange particles for which freeze- 
out may be occurring somewhat earlier. The Tf can be obtained from the intercept in Fig. 
1, while the average collective velocity can be obtained from the slope. 
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Figure 1. r s ] ope for various particles as a function of mass. 

Strangeness changing - chemical- reactions are typically slower than the elastic pro- 
cesses and hence are expected to freeze-out before the thermal freeze-out. The temperature 
and chemical potential at this freeze-out decides the particle yields of various types, pro- 
vided these yields are measured for the full 4ft-integrated region; otherwise the measure- 
ments will depend upon the details of the collective flow mentioned above. Furthermore, 
taking ratios of such yields, one can reduce the dependence on the collective dynamics even 
more. A simple thermal model of free particles at a temperature T, volume V and chemical 
potential p.g has been shown [11] to describe beautifully 22 ratios of particle yields which 
vary by three orders of magnitude, leading to Tf em ~ 170 MeV and p.ff> ~ 270 MeV. Fig. 
2 displays the thermal and chemical freeze-out points for the SPS Pb-Pb collisions along 
with those of other experiments. The figure is taken from Ref. [12] where the references 
for the data are also given. A comment about p.g"J^ al may be in order, as we discussed 
above the corresponding Tf only. Since chemical equilibrium is lost earlier, it is strictly 



4 



Pramana - J. Phys., 



> 

hO 


20 


m 

•3 

*f 
O 


NA4 9 

Nix -Chapman 
Nix et al. 
Dobler et al. 




- <E>/<ll>=lGeV 

- £=45 MeV/fm 3 

- n = 0.05/fm - 










IM 









15 


□ 




*\ 









10 















05 


----- 




to 
i 
to 


S/B= 36 

10 
(5 

«; 


\\ 

0) 
H 

to 











D .2 


0.4 


6 0.8 



,u B [GeV] 

Figure 2. Chemical and thermal freeze-out points in the (T,/jb) plane for various 
experiments. Taken from Ref. [12]. 

speaking not well defined. One simply adjusts f/£ e J" ial such that the particle ratios at Tjrf 7 " 
agree with the observed values. 

Since T^ em turns out to be very close to that expected for the quark-hadron transition 
from lattice QCD, it is plausible that the hadronic chemical equilibrium is a direct conse- 
quence of a pre-existing state of uncorrelated quarks and antiquarks and not due to hadronic 
rescatterings/reactions, since there is not much time for the latter. Hadron formation is then 
governed by the composition of the earlier state in a statistical manner and an expansion 
later does not change their yields. Needless to say though, the proximity of the two tem- 
peratures mentioned above is merely indicative. Indeed such temperatures and chemical 
potentials could still be reached via an expanding hadron gas as well. One then would 
expect though that the particle ratios will not reflect the underlying quark symmetries de- 
scribed next. 



2.3 Excess Strangeness 

Since the early days of heavy ion collisions, when sulphur beams at 6.4 TeV were bom- 
barded on sulphur targets, a global enhancement of strangeness in these collisions has 
been observed relative to e + e~ or pp collisions. Defining [13] a parameter X s to count the 
strangeness, \ s = 2{s + s) / \{u + u) + (d + d)) produced, it was found that A.^ 4 ~ 2Xf /> . For 
both S (6.4 TeV) +Ag and Pb (33 TeV) + Pb, a similar factor of 2 enhancement was ob- 
served. This global enhancement together with the picture of statistical hadron formation 
discussed above suggests an interesting pattern for specific particles. 

Fig. 3 shows the results from WA97 for specific enhancements of CI, S , A etc. The 
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Figure 3. Specific strangeness enhancement in Pb-Pb collisions seen by the WA97 
collaboration [14] as a function of number of participants in the collisions. 



yield for these particles per participant nucleon for Pb-Pb relative to that in p + Be is 
displayed as a function of average number of participants. The number of participants is an 
indicator of how central (or violent) the collision is; maximally central events are expected 
to have largest energy density deposited. One sees the yields to be almost independent 
of the centrality reached experimentally so far. Assuming now that the global strangeness 
enhancement by a factor of two indicates the relative probability of finding a strange quark 
in Pb-Pb vs p — Be to be twice, due perhaps to a formation of a QGP-state in the latter 
which evaporates just before the observed chemical freeze-out, one expects a factor of 
two enhancement for every extra strange quark (or antiquark) in a hadron, explaining the 
pattern in Fig. 3. Note that the masses of these particles increase as more strange quarks are 
added, mn > mn > m\. Their production in a purely hadron rescattering/reactions scenario, 
therefore, will be subject to increasingly higher thresholds, resulting in an opposite pattern 
to that observed. Thus, the pattern in Fig. 3 of specific enhancements of strangeness clearly 
points to quark degrees of freedom. It will be interesting to see whether this enhancement 
sets in smoothly or abruptly as the number of participants vary. 



2.4 Excess Low Mass Dileptons 

An interesting anomaly, observed first in S (6.4 TeV)+ Au collisions by NA45/CERES [15] 
and confirmed for Pb (33 TeV) +Au collisions [16], is an enhancement of M e + e - spectrum 
in the region 250 MeV < M e+e - < 700 MeV. The data for p (450 GeV) +Au collisions in 
the entire range < M e + e - < 1500 MeV can be explained by including the contributions 
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from all known hadron decays of Jt° , r| , co, p , (j), T|' etc. within the acceptance of the detector. 
A similar exercise for S + Au or Pb +Au reveals an enhancement in the low mass region 
mentioned above, with the enhancement factor being 2.6 ± 0.5 (stat) ± 0.6 (syst). Various 
theoretical explanations have been offered [17] to explain it. For us, it is interesting to note 
that thermal emission from an expanding fireball with parameters similar to those in sec. 
2.2. can account for the excess, thus supporting to the picture discussed there. 

2.5 J Suppression 

As remarked in the introduction above, one needs to employ 'harder' probes to explore 
the physics of the fireball at earlier times when QGP may have existed. Production of 
J/\y is one such hard probe. Since it is a tightly bound meson of charm and anticharm 
quarks, Matsui and Satz [18] argued that color Debye screening of these heavy quarks will 
prevent formation of //\|/, if QGP is formed in the heavy ion collisions. Due to a finite 
size and lifetime of the fireball, the observable effect is expected to be a suppression in the 
production of J/\\f. The NA38 and NA50 collaborations [19] measured J /\\f cross sections 
for a variety of collisions, starting from p + d to Pb + Pb using the same muon spectrometer 
in the same kinematic domain (0 < y c J+ M - < 1 and |cos© cs | < 0.5). While the systematic 

errors are thus minimised, the lighter beams were necessarily of high energies; y/Sjw 
thus varies from 17 GeV to 30 GeV. From the observed dimuon spectrum, dN /dM fl + fJ ~, 
obtained after subtracting the background due to like-sign dimuons, the cross section 
was obtained by a fit of the observed spectrum for M^+^- > 2.9 GeV with 5 parameters: 
the number of events of Drell-Yan continuum, //\|/, and \|/' and the J/\\t mass and width. 
Comparing with O^Q th , where isospin corrections were taken into account, the K- 

factor was found to be universal in pp, pA and AB collisions: c® Y B °c A ■ B for all of them, 
where A and B are the mass numbers of the projectile and target respectively. Normalizing 
Bfi + fi~ a AB by dividing by A B therefore, where B^^- is the branching fraction of J/\\t 
in to jU + /j~, one could expect QGP formation to be signalled by a drop at some value of 
A ■ B. Fig. 4 shows the NA38 and NA50 results where one notices a gradual fall in with 
A ■ B for all values. Note that some measurements have been re-scaled so that all are for the 
same energy in this figure. The decreasing cross section for all values of A B, including 
small ones, is an indication of the presence of yet another mechanism for 7/\|/-suppression 
in these collisions. Thus any suppression due to QGP will have to be over and above this 
'normal suppression' . 

Production of heavy quarkonia is an old and mature area of perturbative QCD. In partic- 
ular, hadroproduction of J/\\t has been explained both in the colour evaporation model [20] 
and the colour octet model [21] at y/s comparable to those in Fig. 4. So it is a natural ques- 
tion to ask whether the decrease in Fig. 4 can be explained using pQCD. Unfortunately, 
sufficient information on the nuclear structure functions, especially the gluonic ones, is not 
available at present; assuming them to be independent of mass number A or B is perhaps 
incorrect in view of the famous EMC-effect. Using the existing models of the EMC-effect, 
on the other hand, one finds hardly any decrease in the cross section in Fig. 4. It is even 
likely that this lack of decrease of B^+^-O^g /AB with AB is a generic feature independent 
of the models of the EMC effect. This is due to the fact that the dominant contribution to 
the cross section in Fig. 4 comes from the so-called antishadowing region in x which is 
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Figure 4. J/\\f cross section times its branching fraction in to dimuons as a function of 
product of the mass numbers of target and projectile, AB. The data are from Ref. [19] 
and the la (full lines) and 2c(dashed lines) curves are from Ref. [23]. 

hard to avoid for even the gluons due to the momentum sum rule. In view of the continuous 
decrease in Fig. 4, i.e. even for p+ light-A, where the radius of the target is only 2-4 times 
larger than that for the hadroproduction [20,21], one has to ask whether a pQCD descrip- 
tion of total cross sections for is at all possible. It would be interesting and desirable 
to thrash out this question by extensive investigation of the nuclear glue and its impact on 
the J/\\f cross section. 

The normal suppression in Fig. 4 has been explained [22] as a final state interaction. The 
produced 7/\[/-state or its precursor can get absorbed in the nuclear matter (of the target 
and beam). Treating oj^ as a free parameter and using the known nuclear profiles, one 
finds that a O a b s ~ 6.4 mb can explain the linear fall in Fig. 4 quantitatively in Glauber type 
models. However, the Pb-Pb data point seems to be off this linear fall, and exhibits thus 
an 'anomalous suppression'. One can alternatively use an empirical (AB) a fit to all points 
except the Pb-Pb, which too will be linear on the scales of Fig. 4 and the Pb-Pb data point 
stands out again. 

Unfortunately, the issue of how statistically significant this anomalous suppression is 
gets affected by the crudeness of the theory described above as well as by the assumptions 
needed to rescale some of the data points. Ignoring these systematical theoretical errors, 
one finds the anomalous suppression to be a 5o effect [19], while including them leads 
[23] to a conclusion that no anomalous suppression exist at a 2a or 95% confidence level, 
as shown by the 2o-band (enclosed by dashed lines) in Fig. 4. 

The NA50 collaboration also measures //\[/-suppression as a function of the total pro- 
duced transverse energy Ej. Using the same procedure as outlined above, the number of 
J/\y events and the Drell-Yan events 1 in each £Y-bin are determined. By taking their ra- 



Note that the usual Drell-Yan cross section in pQCD is defined only by integrating all Ej. 
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tio, one obtains a less systematic error prone R expt = B /I + IJ -a J ^ /°m^-m 2 as a function of 
Ej, where M1-M2 is the range of dimuon mass over which the Drell-Yan cross section is 
integrated. Using simple geometrical models, Ej can be related to the impact parameter b 
at which the two nuclei collide. Furthermore, any given b{Ej) can be related to an average 
nuclear path length L which the produced J/\y (or its precursor) has to traverse and which 
will determine the probability of its absorption in nuclear matter. 




2 3 4 5 e "7 a 9 1 o 1 1 

l_(fm) 

Figure 5. The ratio of J/\\l cross section and the Drell-Yan cross section vs. L in fm. 
The crosses are NA38 data, shown along with the straight line fit, a 4a band [24] around 
it, and the NA50 data (squares) with 4a errors on them. The data are from Ref. [19]. 

Fig. 5 shows R expt as a function of L, as determined by the NA50 collaboration, using M\ 
= 2.9 and M2 = 4.5. The normal nuclear suppression can be well approximated by R exvt = 
A ■ exp(— p nuc i ■ O a t,s ■ L) or can be calculated more exactly in a Glauber model. The straight 
line in Fig. 5 displays the fit for the light nuclei for p nuc i —0.17/ fm 3 and c a b s ~ 6.6 mb. 
The low L point for Pb-Pb collisions, corresponding to peripheral collisions, falls on the 
fitted line while all the large L points fall below it. Again, one can ask for the statistical 
significance of this anomalous behaviour. Since the fit above uses data from £V-bins, or 
equivalently L-bins, for lighter nuclei, there are again sizeable errors on the theoretical 
prediction. For the 1995 data, which seem broadly in agreement with the 1996 data and 
the 1998 data, it has been estimated [24] that all the Pb-Pb data points fall in a 4o-band 
although they are all systematically below the theoretical prediction, as shown in Fig. 5. 

It seems thus likely that an additional mechanism to suppress 7/\[/ production in Pb-Pb 
collisions, especially at large values of transverse energy Ef, or large L, is needed over 
and above the normal suppression caused by absorption in the surrounding nuclear matter. 
There have been several theoretical attempts to provide such a mechanism including, of 
course, invoking a possible a quark-hadron transition. A key non-QGP scenario invokes 
the possibility of destruction of the by the so-called co-mover debris of the collisions. 
In a recent 2 such work, it has been claimed [25] that the entire NA50 £y-spectrum of 



2 Although this came after my talk, I include it here for completeness. 
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the J/\\t cross section ratio R can be explained using the co-mover picture. In fact, the 
second shoulder in the Et -behaviour observed [26] in the 1998 data and which could be 
due to QGP [27], has been explained in Ref. [25] as due to fluctuations at the tail of the 
EY-spectrum. The difference between this mechanism and a QGP model will, therefore, 
most likely show up at the upcoming RHIC collider in BNL where Au (19.7 TeV) + Au 
(19.7 TeV) collisions will be studied this year and the Ej tail will extend much farther. 

3. Conclusions and Outlook 

An important non-perturbative prediction of (lattice) QCD is the existence of a new phase 
of matter, Quark-Gluon Plasma, at sufficiently high temperatures. Since the Standard 
Model has so far been tested experimentally only in the weak coupling regime, it seems de- 
sirable to confront this prediction with experiments. Collisions of heavy ions at very high 
energy may be able to deposit the required high energy density over a reasonable volume. 
The experimental programs at BNL, New York and CERN, Geneva have by now provided 
results for Au on Au and Pb on Pb at y/s ~ 0.9 TeV and 3.6 TeV (or y/s NN ~ 5 GeV and 
17 GeV) respectively. The year 2000 should witness a factor of about 39 increase in the 
colliding CMS energy at BNL while LHC at CERN should achieve a y/s = 1 150 TeV The 
experiments so far have provided tantalizing hints of the new phase and therefore of the 
exciting physics in the years ahead. 

A fireball of QGP produced in these collisions cools by expanding and converts into 
ordinary hadrons and leptons fairly quickly. Since this makes a distinction of events with 
QGP formation from those without it a very tough task, it seems prudent to look for a 
congruence of various signatures in as many different ways of detecting QGP as possible. 
Interestingly, the current results do indicate such a trend of congruence from apparently 
unrelated measurements. 

Soft hadron production data, such as /^-distributions, two particle (HBT) correlations 
and ratios of particle yields, can be interpreted in terms of a chemical freeze-out, followed 
by a thermal freeze-out. The freeze-out temperature for the former for the CERN SPS data 
turns out to be <~ 170 MeV ~ T c (quark-hadron transition), suggesting that the hadrons 
were formed from an uncorrected QGP-like state. The global excess of strangeness by 
a factor of two and the specific enhancement pattern of £1, E and A, seen by the WA97 
experiment, showing larger enhancement for the heavier particles with more strange quarks 
is consistent with this picture. In fact, the increasing thresholds make it very difficult to 
explain this pattern in any conventional hadronic picture. Thermal effects arising from T ~ 
120 MeV, the thermal freeze-out temperature, can account for the observed excess of low 
mass dielectron events of the NA45 experiment. 

Finally, anomalous J/\\f suppression seen by the NA50 experiment for Pb-Pb collisions 
can be understood as arising out of a deconfined quark-gluon plasma. Nevertheless, much 
more theoretical and experimental work will be needed to make a convincing case of quark- 
gluon plasma formation in the heavy ion experiments since the signals are still not spec- 
tacular in their statistical significance and credible alternative explanations exist in many 
cases for the observed results . Clearly, the commissioning of RHIC will be a big boost 
and will hopefully result in making a definitive case for quark-gluon plasma. 



10 



Pramana - J. Phys., 



References 

[1] See, e.g., Review of Particle Physics, Eur. Phy. J., C3, 81 (1998). 
[2] K. G. Wilson, Phys. Rev. D10, 2445 (1974). 

[3] M. Creutz, Quarks, gluons and lattices, Cambridge Mongraphs (1983). 

[4] Proceedings of a) Lattice 98, Nucl. Phys. B(PS)73, 1999, b) Lattice 97, Nucl. Phys. B(PS)63, 
1998 etc. 

[5] For a review, see, R. V. Gavai in "Quantum fields on the computer", Ed. M. Cruetz, World 

Scientific 1992, p. 51. 
[6] C. Bernard et al., Phys. Rev. D55, 6861 (1997). 
[7] J. D. Bjorken, Phys. Rev. D27, 140 (1983). 

[8] M. Stephanov, K. Rajagopal, and E. Shuryak, Phys. Rev. D60, 114028 (1999). 
[9] Proceedings of a) Quark Matter 99, Nucl. Phys. A661, 1999, b) Quark Matter 97, Nucl. Phys. 
A638, 1998 etc. 

[10] NA49 Collaboration, T. Alber et al., Phys. Rev. Lett. 75, 3814 (1995). 

[11] P. Braun-Munziger, I. Heppe and J. Stachel, Phys. Lett. B465, 15 (1999). 

[12] J. Cleymans and K. Redlich, Phys. Rev. Lett. 81, 5284 (1998). 

[13] F. Becattini, M. Gazdzicki and J. Sollfrank, Eur. Phys. J. C5, 143 (1998). 

[14] WA97 Collaboration, E. Andersen et al, Phys. Lett. B449, 401 (1999). 

[15] CERES Collaboration, G. Agakichiev et al, Phys. Rev. Lett. 75, 1272 (1995). 

[16] CERES Collaboration, G. Agakichiev et al, Phys. . Lett. B422, 405 (1998); B. Lenkeit in a) of 

Ref. [9]. 
[17] R.Rappina)ofRef. [9]. 

[18] T. Matsui and H. Satz, Phys. Lett. B178, 416 (1986). 

[19] NA50 Collaboration, M. C. Abreu et al, Phys. Lett. B450, 456 (1999); ibid, B410, 337 (1997). 
[20] R. V. Gavai, D. Kharzeev, H. Satz, G. A. Schuler, K. Sridhar and R. Vogt, Int. J. Mod. Phys. 

A10, 3043 (1995). 
[21] S. Gupta and K. Sridhar, Phys. Rev. D54, 5545 (1996). 

[22] A. Capella, J. A. Casado, C. Pajares, A. V. Ramallo and J. Tran Thanh Van, Phys. Lett. B206, 
354 (1988); C. Gerschel and J. Hiifner, Phys. Lett. B207, 253 (1988); Z. Phys. C56, 171 (1992). 
[23] R. V. Gavai and S. Gupta, Phys. Lett. B408, 397 (1997). 
[24] R. V. Gavai, Mod. Phys. Lett. A14, 821 (1999). 
[25] A. Capella, E. G. Ferreiro and A. B. Kaidalov, hep-ph/0002300. 
[26] NA50 Collaboration, M. C. Abreu et al, CERN-EP- 2000-013, Phys. Lett. B, in press. 
[27] S. Gupta and H. Satz, Phys. Lett. B283, 439 (1992); H. Satz, Nucl. Phys. A661, 104c (1999). 



Pramana - J. Phys., 



11 



